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Abstract. We study the distribution of bright star-forming complexes in a homogeneous sample of 72 late-type 
("irregular") dwarf galaxies located within the 10 Mpc volume. Star-forming complexes are identified as bright 
lumps in B-band galaxy images and isolated by means of the unsharp-masking method. For the sample as a 
whole the radial number distribution of bright lumps largely traces the underlying exponential-disk light profiles, 
but peaks at a 10 percent smaller scale length. Moreover, the presence of a tail of star forming regions out to 
at least six optical scale lengths provides evidence against a systematic star formation truncation within that 
galaxy extension. Considering these findings, we apply a scale length-independent concentration index, taking 
into account the implied non-uniform random spread of star formation regions throughout the disk. The number 
profiles frequently manifest a second, minor peak at about two scale lengths. Relying on a two-dimensional 
stochastic self-propagating star formation model, we show these secondary peaks to be consistent with triggered 
star formation; for a few of the brighter galaxies a peculiar peak distribution is observed that is conceivably due to 
the onset of shear provided by differential rotation. On scales between 100 and 1000 pc, and by taking into account 
exponential-disk structure, bright lumps reveal cluster dimensions between 1.3 and 2, with a weak trend to higher 
dimensions for brighter galaxies. Cluster dimension weakly anticorrelates with the lumpiness index (the fraction 
of the total galaxy light due to the light contributed by the lumps), the latter index showing no dependence on 
luminosity. Lump spreading within the disk, as measured by the concentration index, and lump clustering, as 
given by the cluster dimension, are not linked to each other. Interpreting cluster dimension in terms of porosity 
of a self-similar intragalactic medium, we derive a relation between current star formation rate, scale length, and 
porosity. 
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1. Introduction 

Dwarf irregular galaxies exhibit peculiar morphologies 
that are dominated by the flashy but seemingly irreg- 
ular presence of star-forming regions. Unlike with more 
massive disk systems that modulate star formation by 
spiral density waves, dwarf irregulars in the field — or 
non-interacting irregulars in general — constitute ideal 
testbeds for the study of genuine processes regulating lo- 
cal and global star formation and consequently of galactic 
evolution. A review addressing several key questions con- 
cerning large-scale star formation in irregular galaxies is 
given by Hunter (1997), and an evaluation among sim- 
ple models for the onset of star formation in irregulars is 
provided by Hunter, Elmegreen, & Baker (1998). 

Important clues as to hidden constraints shaping the 
heterogeneous appearance of dwarf irregular galaxies may 
emerge from detailed investigations concerning the spa- 
tial distribution of star- forming regions (e.g., Feitzinger 
& Braunsfurth 1984). In recent years several studies on 
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the distribution of H II regions and of young, compact 
star clusters in late-type spirals and in dwarf irregular 
galaxies have appeared (Telles, Melnick, & Terlevich 1997; 
Brosch, Heller, & Almoznino 1998; Elmegreen & Salzer 
1999; Heller et al. 2000; Roye & Hunter 2000; Billett, 
Hunter, & Elmegreen 2002). Applying measures like con- 
centration, asymmetry, and fractional-luminosity indices 
these authors found the star forming regions to be dis- 
tributed rather randomly, with some tendency to central 
concentrations particularly for star-bursting systems. 

In this paper we extend these previous studies by 
analysing the distribution of bright spots or lumps in the 
_B-band images of a sample of 72 late-type ("irregular") 
dwarf galaxies. The general equivalence of bright lumps in 
Ha and in broad band blue images as tracers of star forma- 
tion complexes can be appreciated by comparing galaxy 
images filtered at the two corresponding wavelengths (e.g., 
Elmegreen & Salzer 1999; Sparke & Gallagher 2000, pp. 
139 and 229). With our homogeneous and relatively large 
sample at hand we aim at comparing three morphologi- 
cal indices with each other (as applied to lumps within 
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Table 1. Summary of terms 



Lumps: bright residual features seen in galaxy images after a median filtered version is subtracted from 

the original image; synonymous expressions: bright spots or knots; physical correspondents: star- 
forming complexes encompassing H II regions, young star clusters, and stellar associations; 

Lumpiness index X : fractional flux or ratio of the flux due to the lumps within the residual image and the total galaxy 

light of the original image. 

Cluster dimension D: correlation dimension for a discrete set of lump centers in a plane, i.e. within a radius r around 
a typical lump there are n oc other lumps; no weighting for lump size or luminosity. If consi- 
dered as an indirect measure for a three dimensional medium's fractal dimension, D may be 
related to the volume filling factor of the empty regions, called porosity. 

Concentration index CI: concentration index as the ratio of lump centers in an inner circle and lump centers in an outer 
annulus, normalized according to some prescription; no weighting for lump size or luminosity. 



HTegulars), none of which has been previously reported in 
the form presented here or within our context. Indices of- 
ten serve as the quantitative counterparts to qualitative 
physical concepts. In particular, lump spreading within a 
galaxy may be described by concentration indices of differ- 
ent apertures; lump clustering may be represented by the 
correlation dimension for the two-dimensional lump distri- 
bution; finally, the lumpiness (or flocculency) of a galaxy 
may be measured by means of some fractional light index. 
Table 1 gives a summary of terms and indices that will be 
more carefully introduced in the subsequent sections. 

Comparing relations among morphological indices, we 
may deepen our insights into the various mechanisms re- 
sponsible for the morphology of irregulars. The interstellar 
matter of dwarf irregulars with different global proper- 
ties may be different (e.g., metallicity, mean gas density, 
turbulence, gravitational potential), implying differences 
in the conditions for the formation of stars. Thus aspects 
concerning the abundance and distribution of star-forming 
regions within the galaxies, like clustering and star for- 
mation rates, may turn out to vary correspondingly. We 
adress this issue by means of the clustering parameter. 
Another goal of the paper is to contribute to the discussion 
concerning the influence of shear due to differential rota- 
tion on star formation in gas-rich late-type galaxies (Roye 
& Hunter 2000; Elmcgreen, Palous, Ehlerova 2002). Based 
on cellular automata simulations we introduce a possible 
criterium to be checked for in radial lump number distri- 
butions; we show that a few of our galaxies indeed meet 
the criterium, but more research is needed for conclusive 
results. 

The paper is organised as follows. Section 2 gives an 
overview of the galaxy sample used and presents a table 
with basic galaxy parameters as well as the parameter val- 
ues deduced in the subsequent sections. Section 3 describes 
the adopted lump detection method and introduces a first 
index, the lumpiness index. Section 4 presents the radial 
number and number density distributions for the bright 
lumps of all the galaxies. A concentration index that is 
normalized according to the galaxies being exponential- 
disk systems is introduced and applied to the bright lump 



distribution. We then extensively discuss the possible rea- 
son for a peculiarity seen in the lump number distribu- 
tion, namely the occurence and relative locations of major 
and minor peaks. In Section 5 we determine the cluster or 
correlation dimension for the two-dimensional lump pat- 
terns, relate it to galaxy absolute magnitude, and deduce 
a model mean porosity that is linked to the current star 
formation rate. In Section 6 we check for relations among 
the three indices introduced in the previous sections. We 
end with a discussion and the conclusions in Section 7. 
Image processing was performed throughout within the 
IRAFQ package. 

2. The sample 

The sample of 72 irregular dwarf galaxies is a compila- 
tion of Im, BCD, and late-type S galaxies fainter than 
Mb = — 18 mag — dwarf "irregulars" for short — used 
previously for an analysis of exponential-disk model pa- 
rameters relating photometric, kinematic, and environ- 
mental properties (Parodi, Barazza, & Binggeli 2002). The 
galaxies are lying in the field or in groups within the 
nearby 10 Mpc volume. They were imaged with 1.2-m to 
1.5-m telescopes with resolutions of 0.39-0.77 arcsecs per 
pixel and under seeing conditions varying between 0.8 and 
4.0 arcsecs. Data reduction was performed consistently by 
our team along the usual prescriptions. _B-band galaxy im- 
age galleries can be found in Parodi et al. (2002) and the 
references given therein (or at the bottom of Table 2 in 
the present paper). 

In Table 2 we list basic galaxy data as well as paramet- 
ric data as deduced in the subsequent sections. The first 
nine columns are mostly taken over from Table 1 in Parodi 
et al. (2002), while the seven other columns correspond to 
results obtained in the present paper. The columns read 
as follows: 

Columns 1 and 2 give the galaxy name and the galaxy 

type; 

columns 3 and 4 list the i?-band absolute magnitude and 

^ IRAF is d istributed by the Na tional Optical Astronomy 
Observatories; http:/ /iraf. noao.edu. 
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Table 2. Galaxy and bright-lump data for the 72 irregular dwarf galaxies of our sample 
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(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 
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12 


1.12 
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3.0 


0.37 
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Im 


— 12.66 


22. 15 
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0.50 


22 


I 
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1.0 


0.0 


1.26 
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Im 


— 14.68 


23.89 


1283 
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0.85 


28 


I 


1 












0.019 


UGC4998 


Im 


— 15.78 


22.50 


843 


1879 


0.50 




I 


10 


1.14 






1.13 




0.058 


Hoi 


Im 


— 15.44 


22.75 


1007 


2089 


0.83 


28 


I 


85 


0.90 


2.0 


1.0 


0.81 


1.30 


0.180 


BKIN 


Im 


— 12.70 


24.49 


529 


1832 


0.42 




I 















0.035 


NGC297d 


Sd 


— 17.37 


20.52 


807 


2988 


0.42 


54 


I 


218 


0.36 


2.0 




0.49 


1.59 


0.077 
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— 13.76 


22.30 
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0.67 


27 


I 


4 












0.051 
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— 14.75 


22.73 


593 


1209 


0.56 
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Im 


— 14.18 


24.13 


929 
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1.00 


66 
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0.069 


Kar73 


Im 


— 10.81 


24.53 
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0.67 




I 


1 












0.030 


UGC/8215 


Im 


— 12.71 


22.29 


234 


593 


0.74 


14 


III 

















UJJUlOY 


Im 


— 12.34 


22.64 


248 


565 


0.60 


17 


III 


10 


0.78 


1.5 




0.55 




0.083 


DUUIO8 


Im 


— 14.97 


21.89 


715 


2048 


0.40 


31 


III 


70 


0.58 


1.0 


2.0 


1.47 


1.31 


0.110 


DD0169 


Im 


— 15.30 


22.15 


1032 


2692 


0.32 


26 


III 


15 


0.74 


1.0 




1.92 




0.097 


T T /~i r 

UGC8508 


Im 


— 13.69 


21 .24 


293 


1005 


0.56 


26 


III 


33 


0.96 


2.0 




0.63 


1.82 


0.042 


NGC5229 


Sd 


— 14.44 


20.51 


482 


2007 


0.19 


56 


III 


30 


1.10 


1.0 




1.56 


1.56 


0.120 


ivT /~i r" 

NGC5238 


Sdm 


— 15.03 


21.90 


544 


1507 


0.68 


19 


III 


22 


0.46 


0.5 




6.27 


1.56 


0.086 


DD0181 


Im 


— 13.30 


22.00 


372 


1064 


0.44 


21 


III 


44 


1.11 


0.5 


3.5 


0.99 


1.67 


0.207 


L>1jU18o 


Im 


— 13.90 


21.89 


337 


1016 


0.75 


16 


III 


7 












0.035 


UGC8833 


Im 


— 11.95 


22.41 


177 


425 


0.71 


19 


III 


2 












0.020 


HoIV 


Im 


— 15.95 


22.11 


1634 


4442 


0.26 


39 


III 


65 


1.06 


1.0 


4.0 


0.71 


1.78 


0.175 


1NGG5474 


Scd 


— 17.52 


20.75 


992 


3741 


0.96 


41 


III 


198 


1.41 


1.5 


4.0 


0.50 


1.55 


0.109 


NG05477 


Sm 


— 15.24 


21.59 


559 


1763 


0.70 


32 


III 


37 


1.01 


2.0 


1.0 


0.49 


1.61 


0.121 


DDO190 


Im 


— 15.17 


20.99 


382 


1393 


0.88 


44 


III 


24 


1.02 


1.0 


0.0 


1.23 


1.49 


0.030 


DD0194 


Im 


— 15.00 


23.03 


981 


1814 


0.64 


38 


III 

















UG06541 


Sm/BCD 


— 13.40 


22.49 


294 


677 


0.53 


17 


IV 


7 




0.5 




00 




0.037 


NGC3738 


Irr 


— 15.80 


22.40 


685 


1544 


0.79 


69 


IV 


74 


0.28 


0.0 




12.3 


1.95 


0.062 


NGG3741 


Im/BCD 


— 13.34 


21.71 


207 


609 


0.77 


38 


IV 


14 


1.04 


1.0 


2.0 


1.95 




0.029 


1JJJU99 


Im 


— 14.51 


23.09 


904 


1560 


0.47 


23 


IV 


70 


0.83 


0.5 


2.5 


1.85 


1.69 


0.132 


NGG40o8 


Sm/BCD 


— 15.69 


22.89 


950 


2089 


0.57 


30 


IV 


111 


0.50 


1.5 




0.64 


1.34 


0.243 


INGG41DO 


BCD 


— 14.06 


22.68 


422 


597 


0.66 


17 


IV 


17 


1.14 


0.5 




1.88 




0.076 


UGC7298 


Im 


— 13.71 


22.44 


413 


959 


0.67 


19 


IV 

















JN GG4248 


IBm 


— 16.33 


21.57 


1066 


3313 


0.43 


42 


IV 


16 


1.02 


0.5 




00 




0.072 


DD0127 


Sm 


— 14.39 


24.72 


1003 


1455 


0.51 


32 


IV 


3 












0.039 


UGC7639 


dSO/BCD 


— 15.58 


22.57 


1015 


2184 


0.54 


25 


IV 


64 


0.54 


0.5 




4.24 


1.48 


0.041 


UGG288 


Im 


— 13.82 


22.98 


441 


799 


0.45 


25 


VI 


1 












0.000 


UGG685 


Im/BCD 


— 14.92 


21.96 


509 


1392 


0.63 


38 


VI 


1 












0.019 


UGG1281 


Sd 


— 15.83 


20.90 


1125 


3725 


0.16 


50 


VI 


94 


0.70 


0.5 


2.0 


1.96 


1.69 


0.140 


NGG1156 


IBm 


— 17.68 


20.31 


922 


3971 


0.62 


55 


VI 


63 


0.99 


1.0 


2.0 


0.98 


1.80 


0.118 


UGC2684 


Im 


— 13.13 


23.30 


618 


1100 


0.38 


37 


VI 


1 












0.055 


UGC271d 


Sm 


— 15.08 


23.24 


1047 


1673 


0.47 


27 


VI 


7 




0.0 




00 




0.044 


UGC2905 


Im 


— 14.41 


21.30 


277 


889 


0.63 


26 


VI 


1 












0.016 


UGC3303 


Sd 


— 15.90 


23.07 


1718 


3031 


0.61 


79 


VI 


7 












0.048 


PGC17716 


SBd 


— 16.79 


21.90 


1310 


3207 


0.76 


51 


VI 


5 








4.11 




0.068 


A0554 + 07 


Im 


— 12.25 


23.46 


254 


366 


0.97 


27 


VI 

















UGG3470 


Im 


— 14.26 


21.54 


484 


1477 


0.26 


47 


VI 


5 








3.18 




0.082 


UGG3600 


Im 


— 13.53 


22.88 


672 


1315 


0.29 


39 


VI 


2 












0.049 


NGC2337 


IBm/BGU 


— 16.39 


21.59 


801 


2637 


0.62 


74 


VI 


33 


0.91 


1.0 


2.0 


2.15 


1.66 


0.094 


UGC4115 


Im 


— 13.51 


22.66 


493 


1054 


0.44 


44 


VI 


16 


0.84 


1.5 


3.5 


0.00 




0.077 


NGC2537 


Sm/BCD 


— 16.60 


22.10 


711 


1911 


0.91 


102 


VI 


53 


0.47 


1.0 




3.34 


1.59 


0.079 


DD064 


Im 


— 13.95 


22.56 


636 


1532 


0.39 


42 


VI 


17 


1.83 


1.5 




0.65 




0.130 


T71 / \ /I "vo r\ ^ 

hjbU47o — G()24 


Im 


— 13.74 


22.00 


454 


1 179 


0.46 


17 


VII 


22 


0.34 


1.5 




0.10 


1.45 


0.196 


ES0115 — G021 


Sm 


— 15.18 


21.63 


1 164 


3455 


0.15 


50 


VII 


372 


0.85 


1.0 




1.24 


1.55 


0.257 


ES0154 — G023 


Sm 


— 16.23 


20.10 


1141 


5635 


0.22 


56 


VII 


318 


1.54 


2.0 


3.5 


0.39 


1.62 


0.204 


IVT /'"^ 1011 

JNGGloll 


Sm 


— 15.69 


21.74 


884 


2587 


0.30 


42 


VII 


29 


0.56 


0.5 




4.02 


1.78 


0.107 


IG1959 


Sdm 


— 15.92 


19.25 


563 


2485 


0.24 


58 


VII 


225 


2.29 


2.5 


0.5 


0.28 


1.56 


0.142 


IC2038 


Sd 


— 14.47 


21.63 


596 


1931 


0.29 


33 


VII 


33 


0.41 


0.5 




4. 10 


1.76 


0.066 




Siii/BCD 


— 16.25 


21.64 


683 


2086 


0.57 


39 


VII 


59 


0.55 


0.5 




7.78 


1.60 


0.083 


AM0521 — 343 


Im 


— 14.24 


21.60 


358 


1020 


0.67 


42 


VII 


11 


1.14 


2.0 


1.0 


0.49 




0.058 


ES0555 — G028 


Im 


— 13.79 


23.43 


870 


1399 


0.37 


54 


VII 


7 




1.0 




1.11 




0.065 


ES0489 — G056 


Im 


— 12.30 


23.00 


186 


377 


0.72 


19 


VII 


22 


0.78 


1.0 




0.82 


1 .47 


0.094 


ESO490 — G017 


Im 


— 14. 18 


21.30 


269 


918 


0.79 


28 


VII 


36 


0.90 


0.5 


1.5 


4.21 


1 .47 


0.076 


ESO308 - G022 


Im 


-13.71 


23.81 


609 


683 


0.90 


33 


VII 


11 


1.26 


0.5 


2.0 


0.46 




0.117 


PGC20125 


Im 


-13.46 


23.75 


676 


502 


0.70 


36 


VII 


61 


0.96 


1.5 


0.0 


2.16 


1.29 


0.292 


ES0558 - PNOll 


Im 


-16.30 


22.08 


825 


2339 


0.71 


68 


VII 


49 


0.25 


1.0 




18.7 


1.38 


0.120 


ESO059 - GOOl 


Im 


-14.49 


22.20 


536 


1288 


0.71 


46 


VII 


60 


0.70 


1.0 




1.20 


1.70 


0.149 


ESO006 - GOOl 


Im 


-14.93 


22.12 


504 


1271 


0.85 




VII 


22 


0.41 


0.0 




9.97 


1.36 


0.104 


UGCA148 


Im 


-14.09 


21.40 


342 


1158 


0.72 


36 


VII 


6 




1.5 




0.34 




0.038 


UGCA153 


Sm/BCD 


-14.21 


23.32 


937 


1427 


0.42 


52 


VII 


23 


0.50 


0.5 




0.77 


1.44 


0.187 


NGC2915 


Sm 


-16.61 


21.26 


673 


2290 


0.56 


62 


VII 


81 


0.36 


0.0 




7.03 


1.98 


0.091 


UGCA193 


Sdm 


-15.15 


21.98 


1480 


4138 


0.11 


55 


VII 


23 


0.56 


0.5 




2.36 


1.24 


0.284 



References to the published papers of our series on Structure aiid stellar conLe7ii of dwarf galaxies: I - Bremnes et al. (1998); /// - Bremnes 
et al. (1999); IV ~ Bremnes et al. (2000); VI - Barazza et al. (2001); VII - Parodi et al. (2002). 
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Fig. 1. ESO 473-G024. Left: B-band image, taken at the 1.5-m Danish Telescope at La Silla, Chile. Right: Residual 
image after subtracting from the original image its median filtered version. 



— from fitting an exponential law to the observed sur- 
face brightness profile — the extrapolated central surface 
brightness, both corrected for galactic extinction; 
columns 5, 6, and 7 give the i?-band disk scale length 
along the semi-major axis (in parsecs), the radius of the 
25th-mag/arcsec^ isophote (in parsecs), and the axis ratio 
of the elliptical isophote, respectively; 
column 8 lists the rotational velocities of the galax- 
ies as far as measured. The four entries for the south- 
ern hemisphere galaxies ESO473-G024, IC2038, ESO490- 
G017, and ESO059-G001 are based on HIPASS | public 
data, treated as in Parodi et al. (2002); 
column 9 points to the five published papers of our series 
with original and detailed observational and photometric 
data on the galaxies. Cf. the reference notes on the bot- 
tom of Table 2; 

column 10 indicates the number of bright lumps found in 
residual galaxy images (Section 3); 

columns 11, 12, and 13 give the scale length of the radial 
number density distribution of the bright lumps (Section 
4.2), the radius of the highest peak in the radial lump 
number distribution, and the peak radius of a possible sec- 
ondary, minor peak, respectively (Section 4.4) (all lengths 
in units of scale length i?^); 

columns 14, 15, and 16 list the values for the following 
indices: the normalized concentration index CI measured 
using i?25 as the outer aperture radius (Section 4.3), the 
reduced cluster (or correlation) dimension D (Section 5), 
and the lumpiness index X (Sections 3 and 6), respectively. 



2 The HI Parkes All-Sky Survey, or HIPASS, is a 21-cm 
HI survey of the southern sky undertaken with a multi- 
beam receiver on the Parkes telescope in Australia. The one- 
dimensional spectral data for a freely choosable position is 
available for downloading in a variety of different form ats at 
http:/ /www. atnf.csiro.au/research/multibeam/release/. 



3. Lump detection and lumpiness index 

A residual image highlighting bright star-forming com- 
plexes — comprising HII regions and young clusters — 
was constructed by subtracting a smoothed version from 
the original B-band image for each galaxy. To obtain the 
former image a convenient method is to median filter the 
original image with a sliding square window. Relying in- 
stead on adaptive filtering techniques (i.e., on the IRAF 
task ADAPTIVE) did not improve on the results, thus we 
kept applying the simple median filtering method. For the 
window w some characteristic metric size should be cho- 
sen: we adopted w — 0.2i?cff, with Rcs being the effective 
radius of the galaxy. An example of images processed this 
way can be seen in Fig. |l|. The galaxy shown is ESO 473- 
G024, with the left frame containing the original image, 
and the right frame presenting the residual image. From 
these latter, high-spatial frequency images we extract the 
following information: 

(i) we detect and tabulate the locations of the bright 
lumps, thus constructing a data base for the distribution 
analysis below. Detection of bright lumps was done au- 
tomatically and thus consistently for all galaxies. First, 
the sky-subtracted residual image is cleaned around and, 
concerning obvious foreground stars, within the galaxy. 
The few foreground stars that probably went undetected 
should have no potential to fake the statistical outcomes of 
this paper. Second, an appropriate point spread function 
(PSF) was looked for: a routine checks for the maximum 
pixel value corresponding to the brightest lump and deter- 
mines its PSF. It may happen that the brightest lump 
takes part in a compact cluster of lumps; thus if the PSF's 
fuU-width-at-half-maximum was more than 10 pixels the 
routine assumed blending and searched for a smaller lump 
with a narrower PSF. Similarly, if the PSF was very peaky, 
an overseen foreground star was assumed and a broader 
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PSF was applied for the lump search. Third, this PSF 
is then used by the IRAF task DAOFIND to search for 
all other lumps above a detection threshold of 3 sky sig- 
mas. In ESO 473-G024, for example, 22 bright lumps were 
thus counted. Some galaxies are too fuzzy-looking to have 
any lumps detected. Finally, lump coordinates are stored 
in physical units with the center and the major axes of 
the 25th-mag/arcsec^ isophotal ellipse providing the ori- 
gin and the axes of the coordinate system; assuming ax- 
isymmetry, the lump coordinates will be deprojected to zero 
ellipticity for all applications. Note that we do not esti- 
mate lump luminosities in order to compare lump bright- 
nesses or to provide luminosity functions; for blue lumi- 
nosity functions of star-forming complexes in spiral and 
irregular galaxies see Elmegreen & Salzer (1999). 

(ii) We determine the galaxy's lumpiness index X ■ 
Lumpiness (or flocculency) of a galaxy may be quanti- 
fied by the high spatial frequency power X , pioneered 
by Isserstedt & Schindler (1986) and recently applied by 
Elmegreen & Salzer (1999) to spiral and irregular galaxies 
and by Takamiya (1999) to HDF-N galaxies with the aim 
of having a galaxy structural parameter at hand that is 
related to the current star formation rate. Following the 
notation of Takamiya (1999) it is defined as the ratio of 
the flux from the bright lumps, and the total flux fg of 
the galaxy, thus 



By measureing the total fluxes of the residual and the 
original images for each galaxy, we obtained the X val- 
ues listed in Table 2. The lumpiness index will be applied 
below in Section 6. 

4. Radial distributions of detected bright lumps 

4.1. Number distribution 

In Fig. ^ we show for all galaxies of the sample the binned 
radial number distributions of the detected bright lumps. 
Histogram bins correspond to concentric elliptical annuli 
with semi-major axes successively growing by half a scale 
length Rd and shown out to 8 Rd', the last bin included 
in the panels comprises all the detected lumps at radii 
larger than 7.5 Rd- The number of counts per bin, A'norm, 
is normalized by the largest bin value within 7.5 Rd', the 
largest bin value as well as the total number of counts are 
printed within each panel. This kind of normalization was 
imposed to minimize the effects of different resolutions 
and seeing conditions for the following inquiry. 

The accumulated radial number distribution N^°l^^ 
for the lumps of all the galaxies, i.e. summing up all 
the profiles of Fig. |^, is shown as histogram in Fig. ^. 
The bright-lump distribution of many galaxies thrown to- 
gether is represented by a radius-weighted exponential 
distribution that is indicated by the solid line obeying 
N{R) oc Rexp{-R/Ri) with Ri = 0.86i?d (cf. Section 
4.2). This basically reflects the exponential light profiles of 



dwarf irregular galaxies in general with, however, a slightly 
shorter scale length than is seen for the -B-band continuum 
light (dashed line) . Note that bright lumps or star- forming 
complexes not only are found way out to large radii, but 
that they constitute a nice tail in the radial number dis- 
tribution out to at least six optical scale lengths, a point 
we come back to in the discussions of Section 7. 

The similar exponential structure for two components 
of the disk is comparable to H II region distributions in 
other types of exponential-disk galaxies. In intermediate- 
type spirals (Athanassoula et al. 1993) and in irregular 
galaxies (Hunter et al. 1998) the azimuthally averaged ra- 
dial distribution of H II regions follows the stellar light 
distribution as well. In our dwarf irregular galaxies bright 
lumps are exhibiting this same behaviour, indicating them 
to be representative for the distribution of H II regions, as 
expected. We can put this statement on a still firmer basis 
tracking down also the radial number density distribution 
of the bright lumps. 

4.2. Number density distribution 

Radial number density distributions for the bright lumps 
are obtained by dividing the number of counts in a given 
bin by the surface of the corresponding elliptical annulus 
(cf. Hodge 1969, Athanassoula et al. 1993). The bright- 
lump number density profiles for our galaxies are shown 
in Fig. ||. Solid lines represent exponential fits for galax- 
ies with a total of at least 10 detected lumps; their scale 
lengths Ri are given in each panel in units of _B-band con- 
tinuum light scale lengths Rd and are listed in Table 2. 
For the whole sample we find a mean of 

< >= 0.86 ±0.06, (7 = 0.41. (1) 

Thus, while the number density distributions of bright 
lumps are — at least partially — also exponential, their 
slopes are on average steeper than those for the density 
distributions of the underlying stellar light. Wc note that 
this very same quantitative behaviour was also observed 
for the number density distribution of H II regions in 
intermediate-type spiral galaxies: Athanassoula et al. 
(1993) found < Rh ii/Rd 0.8, a = 0.4. 

Concluding this subsection we state that dwarf irreg- 
ular galaxies show azimuthally summed-up bright-lump 
profiles that are quantitatively comparable to those of H II 
regions in exponential-disk systems. Thus, as expected, 
bright star-forming complexes largely represent H II re- 
gions. In particular, the scale lengths for H II regions, for 
bright lumps on i?-band images, and for B- and i?-band 
continuum light images (Parodi et al. 2002) obey on aver- 
age the ratio equation Rhii ■ Ri ■ Rd{B) : RdiR) ~ 0.8 : 
0.9 : 1.0 : 1.1, i.e. the older the underlying population 
the larger the scale lengths. This general trend in star- 
forming dwarfs has been observed before, and we discuss 
some implications in Section 7. 
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Fig. 2. Radial number distribution of bright lumps in 72 dwarf irregular galaxies, based on B-band images. Bin width 
is 0.5 Rd,; the rightmost bin sums all counts beyond 7.5 Rd- For each galaxy the counts are normalized by the highest 
bin value. In each panel the galaxy name, the highest bin value, and the total number of counts are given. 



4.3. Concentration index 

The concentration index CI of a galaxy is a convenient 
parameter to quantify galaxy morphology of low- and 
high-redshift galaxies. It compares the light content for 
different radial intervals. Various definitions have been 
used, none of which takes into account the exponential- 
disk constraint. It is, however, a trivial observation that 
in scale length-versus-luminosity diagrams the scatter in 



scale length around the mean relation correlates with the 
concentration index. To have a scale length-independent 
concentration index we will explicitely factor in this un- 
derlying disk feature. Following Heller et al. (2000) for 
the sake of comparision, the concentration index C'I{R) is 
taken as the ratio of the flux or, in our case, the number of 
complexes within an elliptical aperture of semi-major ra- 
dius R/2, i.e. from the inner part of the galaxy, to the flux 
or number of complexes from its outer annulus with inner 
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Fig. 3. Radial number density distribution of bright lumps in 72 dwarf irregular galaxies, deduced from the number 
distribution shown in Fig. |[ In each panel the galaxy name and the distribution's approximate exponential-fit scale 
length are given (in units of B-band scale lengths; evaluated only if there is a total of at least 10 detected lumps, and 
represented by straight solid lines). 



and outer semi-major radii of i?/2 and R, respectively. 
Opposite to Heller et al. (2000) — who suspect linear ra- 
dial distributions of the (flux from) star forming regions 
— we do not bring the two numbers to an equal-area ba- 
sis by dividing the outer number of lumps by a factor of 
three. Instead, we want to relate the measured lump CI 
to the corresponding one for an assumed underlying radial 
number distribution N{R) (x Rexp(—R/ Rd) (or, as it is 



more adequate for lumps, oc Rexp{—R/ Ri); see below). 
We thus normalize our galaxy concentration index with 
C/o(i?) = J^^^ N{r)dr/ j"^^N{r)dr. Expressing the to- 
tal radius R in terms of the disk scale length Rd, R = x Rd, 
one finds 



Cloix) 



[1- (l + a;/2)e-"/2] 
[(1 -hx/2)e-^/2 _ (i + x)e- 



(2) 
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Fig. 4. Total radial number distribution of bright lumps 
in 72 dwarf irregular galaxies. Binning in units of a tenth 
of a scale length; the contributions from each galaxy are 
normalized by the value of its highest-value bin (cf. Fig. 
The dashed line follows a distribution Rex-p{—R/ Rij) as 
it is expected from the continuum light exponential disks. 
A better match is provided for a distribution with a mean 
lump scale length of Ri — 0.86Rd, plotted as solid line. 



Thus our normalization constant CIo{x) depends on the 
chosen outer radius R, or x, respectively, used to deter- 
mine the concentration index. On a statistical footing the 
72 dwarf irregular galaxies of our sample share the mean 
ratios Rd ■ Reff ■ R25 '■ -R995 ~ 1 : 1.5 : 3 : 5.1 in the con- 
tinuum S-band light, where Reff is the effective radius, 
i?25 is the 25th-mag/arcsec^ isophotal radius, and i?995 is 
the radius equivalent to an aperture with 99.5% of the to- 
tal flux. Theoretically, pure exponential-light decays yield 
Rd ■■ Reff ■ R25 ■ R995 = 1 : 1.7 : 4.1 : 7.4, where the 
value for R25 was interpolated using the observed mean 
ratios given above. The fact that our galaxies exhibit on 
average smaller radii is a consequence of the exponential- 
disk description being only an approximation; in particu- 
lar, there is a large fraction of bright galaxies with cen- 
tral light cusps. The appropriate normalization constants 
corresponding to the above radii are then C/o(l)=0.52, 
C/o(1.5)=0.64, C/o(3)=1.23, and C/o(5.1)=3.01. 

With this normalization we expect the concentration 
indices to be CI ~ 1 . Deviations from this canonical value 
provide information on the actual shapes of the profiles. 
For example, the median values of the concentration in- 
dices obtained by Heller et al. (2000) for their dwarf ir- 
regular (BCD + LSB) galaxy sample are — corrected 
to our normalization by adopting their limit R = R25 
and thus applying CIq = 1.23 — (4.23/3)/1.23 = 1.15 
and (3.43/3)/l. 23=0.93 for the Ha-flux and the contin- 
uum images, respectively. Thus BCD and LSB galaxies 
are rather well represented by exponential light profiles on 
large parts. However, the values for the BCDs typically lie 
above and those for the LSBs below these median values, 
reflecting the fact that the former galaxies are more ac- 
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Fig. 5. Normalized bright-lump concentration indices CI 
vs. galaxy rotation velocity. Pure exponential-disk distri- 
butions correspond to CI = 1. Filled circles represent 
galaxies of our sample, triangles at panel upper bound- 
aries stand for infinite CI values resulting from galaxies 
with no lumps detected in the outer annulus, and open 
symbols is data from Roye & Hunter (2000). 



tively star forming in the center regions than the latter. 
Returning to our sample, we arrive at values CI =0.94, 
1.26, 1.32, 1.73 for outer radii R —2, 3, 4, 5 Rd, respec- 
tively. These values typically being larger than one and 
even increasing with larger outer radii is due to the nor- 
malization used so far that was based on the continuum 
light scale length Rd- 

However, if the lump scale length Ri is used for the 
normalization instead of Rd, one indeed recovers CI ~ I. 
This is shown in Fig. 5 where we have plotted the normal- 
ized concentration indices CI{x) of our galaxies for various 
aperture radii R — xRi against their rotational velocities, 
now adopting the mean scale length Ri — 0.86Rd found 
for the lumps in Section 4.2. Actually, one would prefer 
to adopt for each galaxy its particular lump scale length, 
but given the uncertainties in determining them, we are 
content with the mean value given in equation (1). Only 
galaxies with at least five detected lumps are included, 
leaving about 50 galaxies. Infinite values result for CIs in 
the case of no outer-annuli lump detections; in the plots 
they are included as triangle symbols with values fixed 
at CI — 13. Data for the upper right panel is listed in 
Table 2. The panel's median CI values, plotted as dashed 
lines and ignoring the CI=oo cases, are 0.83, 1.24, 0.95, 
and 1.04. Lying all in the vicinity of one, this is consistent 
with the annulus-integrated exponential distribution for 
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the summed lump number distribution seen in Fig. 4 out 
to large radii. 

Thus, applying the analytic tool of the concentration 
index, we have again demonstrated that the radial dis- 
tribution of star forming regions is non-linear but follows 
an annulus-integrated exponential distribution. This im- 
plies a non-uniform random spread of the star forming re- 
gions throughout the disk, which explains the discrepancy 
found by Heller et al. (2000) between the Clna values for 
actual galaxies and the lower ones for simulated galax- 
ies with random star formation region positions. Roye & 
Hunter (2000) pointed out an increased scatter of concen- 
tration indices for faster rotating galaxies of their sample. 
In the upper right panel, data from Roye & Hunter (2000), 
adopted to our normahzation using C/o(3.5) = 1.53 (as- 
suming R25—3Rd=S.5Ri) and with a median CI value of 
only 0.58, are plotted, too; however, we no longer see this 
effect in any panel with our larger sample. 

4.4. Peak number distribution: shear-enhanced star 
formation at work? 

For the sample as a whole the radial locations of bright 
lumps mirror the exponential intensity distribution of the 
underlying population. This is, however, only a rule-of- 
thumb. Individual galaxies may exhibit strong deviations 
from this mean statistical behaviour (cf. Fig. ||). For ex- 
ample, the primary peak in the radial number distribu- 
tion (i.e. the bin with the highest value) appears not 
around one scale length, but is shifted to smaller or higher 
radii. The histogram for the radial distributions of primary 
peaks only is shown in Fig. ^, left (thin-lined bars). Only 
galaxies with peaks corresponding to at least three counts 
are included. The expected maximum of occurences of pri- 
mary peaks around one scale length is clearly recovered; 
a Gaussian with a mean at one scale length and a stan- 
dard deviation of 1.4 is overplotted as the dashed line. 
However, another particular feature of individual lump 
number histograms is the frequent presence of a secondary 
peak that is lower than the main peak (instead of mono- 
tonicly smaller bin heights on both sides of the main peak); 
this is the case for about a third of our galaxies (cf. Fig. 
§). The radii of primary peaks (with at least three counts) 
and secondary peaks are listed in Table 2; in a few cases of 
equal height peaks we refered to 0.1 i?d-bin width number 
distributions to decide which of the peaks is the primary or 
the secondary one. While the primary peaks of those galax- 
ies exhibiting a second, minor peak as well are crowded 
around one scale length (thick-lined bars), the correspond- 
ing distribution of secondary peaks reveals a pronounced 
maximum at about two scale lengths (shaded bars). 

Is this excess of bright stellar complexes at radii larger 
than about two scale lengths a statistical fluctuation or 
is it a manifestation of some underlying mechanism? A 
candidate mechanism that deserves closer inspection is 
shearing due to differential rotation within the outer part 
of the disk. We thus discuss the physical plausibility for 



the influence of shearing on the generation of star-forming 
complexes within dwarf galaxies. Two questions will be 
adressed: First, is shearing in dwarf galaxies a viable 
mechanism? And second, can it account for the observed 
peak distribution? 

A compilation of 20 high quality dwarf galaxy rotation 
curves by Swaters (2001) shows them to look much like 
those of spiral galaxies, with rotation curves rising steeply 
in the inner parts and flattening in the outer parts. In par- 
ticular, most dwarf galaxy rotation curves start to flatten 
around two disk scale lengths, and no dwarf galaxy shows 
solid-body rotation beyond three disk scale lengths any- 
more. Concerning our observed occurence of minor peaks 
in the bright lump distribution, starting at and being 
most pronounced at about two scale lengths as well, we 
may wonder whether it is related to the transition from 
solid-body to differential rotation. Affirmative signals ar- 
rive both from theory and simulations, (i) Larson (1983) 
suggested that the SF rate increases with higher shear 
rate via the "swing amplifier" mechanism: citing Toomre 
he points out that shear itself contributes strongly to the 
growth of gravitational instabilities, leading to gas den- 
sity enhancements and subsequent star formation. For dis, 
however, lacking spiral-density waves, swing amplification 
may seem an inappropriate mechanism to rely on. (ii) 
Alternatively, in their review Seiden & Schulman (1990, p. 
40) remind that in models for stochastic self-propagating 
star formation (SSPSF) shearing increases the density of 
star-forming regions: gas-rich, potential star-forming re- 
gions are transported to and mixed with former star- 
forming regions, giving space for new star formation, (iii) 
Additionally, while it has been questioned whether shear 
may cause any visible effect at all given dwarf galaxies be- 
ing rather slowly rotating systems (Hunter et al. 1998), it 
becomes more and more evident that some irregular galax- 
ies like the Large Magellanic Cloud or NGC4449 possess 
regular, large-scale magnetic fields (e.g., Otmianowska- 
Mazur et al. 2000), and thus it similarly becomes feasi- 
ble that the magneto-rotational instability (e.g., Balbus 
& Hawley 1998) comes into play, effectively strengthen- 
ing the effects of shear. Bearing in mind this possibility, 
we nevertheless restrain the discussion in the following on 
the second of these scenarios only. 

Self-propagating star formation is observed with many 
galactic as well within many extragalactic objects. It is 
thought of as a locally important SF triggering mechanism 
in all types of galaxies. For example, modulated by density 
waves, long-lived spiral arms may be formed in bright disk 
galaxies (Smith, Elmegreen, & Elmegreen 1984); the sur- 
face filling factor of bubbles and the locations of molecular 
rings in observed disk galaxies can be quantitatively ex- 
plained by SSPSF (Palous, Tenorio-Tagle, & Franco 1994); 
age gradients in star-bursting galaxies can be accounted 
for by means of triggered star formation (Thuan, Izotov, & 
Foltz 1999; Harris & Zaritsky 1999); and last but not least, 
the general burst characteristics of compact and irregular 
dwarf galaxies is long known to partially be understood 
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Fig. 6. Left: Peak number distribution of bright lumps in dwarf irregular galaxies, deduced from the number distri- 
butions shown in Fig. H; included are only galaxies with peaks corresponding to at least three counts. The histogram 
for all the primary (i.e. highest) peaks is shown as bars bordered with thin lines; the dashed line is a Gaussian eye fit 
forced to peak at one scale length. The histogram for the primary peaks of those galaxies exhibiting secondary peaks 
as well is overplotted with thick lines; the shaded bars indicate the distribution of the corresponding secondary peaks. 
Note the pronounced signals at the center and in particular at two scale lengths. Right: Peak number distribution for 
20 simulated galaxies, once for galaxies with solid-body rotation (upper panel) and once for galaxies with a transition 
to strong differential rotation at two scale lengths (lower panel). The simulated galaxies are generated by means of a 
stochastic self-propagating star formation (SSPSF) disk model. As in the figure to the left, thick-lined bars represent 
primary peaks whereas shaded bars inform on secondary peaks. Note the reversed peak distribution. 



by means of SSPSF (Gerola, Seiden, & Schulman 1980; 
Comins 1984). 

Relying ourselves on a two-dimensional SSPSF model, 
we numerically tested the hypothesis that the onset of 
shear-induced star formation around the turnover radius 
may leave its imprint in an overabundance of SF regions 
or of bright stellar complexes beyond two scale lengths. In 
the Appendix we describe the particular implementation. 
A general finding of our simulations is that the inclusion 
of shear (i.e. rotation) allows for about five to ten per- 
cent more star forming cells, the exact value depending 
on the particular parameters used. Being mainly inter- 
ested in the azimuthally summed-up radial distribution of 
lumps under different rotational conditions, we compare 
simulation runs with and without a transition to a flat 
rotation curve. In the top panel of Fig. |^, right, a typi- 
cal outcome for a simulation of 20 galaxies with rigidly 
rotating disks is plotted. The highest peaks are found to 
be located around one scale length, while the secondary 
peaks show occurences at many radii but with a preference 
for locations around two scale lengths. This coincides with 
our observed peak distributions. For comparision, in the 
bottom panel of Fig. |^, right, a representative peak dis- 
tribution for a simulation run of 20 galaxies that exhibit 
a continuous transition from solid-body to (strong) differ- 
ential rotation at two scale lengths is shown. Interestingly, 
the primary peaks now occur preferentially at around two 
scale lengths indicating a strong influence of shear on star 
formation around the turnover radius. While this is not 
the general picture observed with our sample, some of the 



brighter galaxies actually do match this pattern: IC 1959, 
ESQ 154-G023, Ho I, NGC 5477. 

We thus conclude that the observed pattern of pri- 
mary peaks at one scale length manifests the underlying 
exponential-disk systems, and that the frequent occurence 
of secondary peaks at about two scale lengths is not neces- 
sarily related to the onset of strong shear in rotating disks. 
As the simulations show, it is however consistent with the 
idea of triggered star formation based on a stochastic self- 
regulation scenario. Some of the larger galaxies are ex- 
hibiting pronounced primary peaks at two scale lengths 
but show minor peaks at one scale length; with these 
galaxies wc may be directly witnessing shear-induced star 
formation. The possible role of bars will be reflected in 
Section 7. 

5. Clustering properties of bright lumps: cluster 
dimensions on scales of a few 100 pc 

Star-forming complexes in dwarf galaxies form non- 
random point patterns also in a sense different from that 
discussed in Section 4.3. Their positions correlate accord- 
ing to a self-similar (fractal) arrangement. In this Section 
we substantiate this claim studying an index devoted to 
spatial statistics, namely the correlation or clustering di- 
mension, as applied to two-dimensional bright-lump dis- 
tributions. 

As observed by Elmegreen & Elmegreen (2001), the 
distribution of bright-lump center positions on a kilopar- 
sec scale in spiral and irregular galaxies obey a power-law 
behaviour similar to the fractal structure of the interstel- 




Fig. 7. Left: Cumulative number vs. aperture radius for the detected lumps within selected galaxies. Logarithms are 
given to base v^. The slope of the straight parts gives the cluster dimension. Fat lines represent large-scale length 
(and high-lump number) galaxies, thin lines short-scale length (and low-lump number) galaxies. Right: Reduced 
bright-lump cluster dimension vs. extrapolated central surface brightness. The symbol size indicates the number of 
lumps used for the determination of the observed (non-reduced) cluster dimension. The line corresponds to a bisector 
fit to the data (equation given in the text). 



lar gas with fractal dimension D3 = 2.3. Thus the center 
positions of star-forming aggregates within isolated areas 
of large galaxies are fractal. Here we address the question 
whether star-forming complexes that are scattered over 
the entire disks of dwarf galaxies are non-randomly dis- 
tributed as well. We restrict our inquiry to the dwarfs of 
our sample that exhibit more than 20 bright lumps. Given 
our photometry with image scales of typically well above 
10 parsecs/pixel and seeing conditions of a few pixels we 
expect to only dissolve structures larger than about 100 
parsecs. Thus small-scale clustering and the accompagn- 
ing blending effects (Elmegreen & Elmegreen 2001) are of 
no concern to our study. To quantify the spatial clustering 
of the position patterns we adopt the cumulative distance 
method (Hastings & Sugihara 1996): a power law relation- 
ship N{r) cx r^ is assumed for the cumulative number of 
points N{r) within a distance r around each point. If the 
distribution is (at least partially) self-similar this will be 
manifested in a log{N)-log{r) diagram as a straight line 
with slope D, called the cluster (or correlation) dimension, 

d log N{r) 
d logr 

The more highly clustered the points (at all relevant 
scales), the lower the cluster dimension. For a random 
or Poissonian distribution of points on a two-dimensional 
plane one has f « 2, independent of the number of 
points involved, which only governs the error estimate. 
The graphs for six observed galaxies with 20-30 lumps and 
for five galaxies with about 200-300 lumps are shown in 
Fig. 7, left, plotted with thin and thick lines, respectively. 
For both groups the relevant scaling range, i.e. the straight 
part of the curve, lies between about 100 (« a/2 ) and 



1000 (« V2 ) parsecs. The galaxies with lower lump num- 
bers exhibit smaller cluster dimensions (£) « 1.5) than the 
galaxies with many detected lumps {D « 1.9). However, 
plotting D versus Niumps for all our data (not shown), 
no clear relation between the two quantities is seen any- 
more. There nevertheless is a hidden dependence between 
the two variables: it emerges from the non- uniform dis- 
tribution of lumps in exponential-disk systems (as dis- 
cussed in Section 4), and it is to be corrected for. We 
do so by, first, simulating point patterns with exponential 
radial number density distributions and indeed are recov- 
ering the observed dependence of the cluster dimension 
on the number of lumps. In particular, accepting a lin- 
ear regression we obtain Dgimui = 0.0013Af;„mps + 1-471. 
Actually, a function converging asymtotically toward D=2 
for large lump numbers would be more appropriate; hav- 
ing no clue as to its exact form, though, we stay within the 
linear approximation. Then, second, instead of using the 
observed chistcr dimensions as inferred from galaxy im- 
ages, we introduce reduced cluster dimensions defined by 
D = Dobs — 0-00^3Niumps, i-e. all measured cluster dimen- 
sions are made comparable by formally adjusting them to 
the common number Niumps=0- Other values could have 
been chosen; however, the adopted value (or other low 
values, say A^i«mps~50) yields consistently cluster dimen- 
sions of about or below the theoretical maximum value 
of two. We show in Fig. 7, right, the reduced cluster di- 
mension as a function of the extrapolated central surface 
brightness for all galaxies. There is a weak but signifi- 
cant trend that fainter dwarf galaxies exhibit lower cluster 
dimensions, i.e. more strongly clustered star-forming re- 
gions, than brighter dwarf galaxies. The same statement 
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holds if instead of central surface brightness we take the 
absolute magnitude of the galaxy. 

Wc have also determined the cluster dimensions for 
15 selected sub-galactic areas (consisting of about 30 
lumps within a circle of about 1.5 kpc diameter) within 
larger galaxies (each with a total of more than about 200 
lumps). With a mean of and a scatter of only 

about 0.1, these areas show relatively high cluster dimen- 
sions that are typically lying above their galaxies' values. 
It furthermore implies that cluster dimensions for local 
lump aggregates scatter less than those for entire galaxies. 

We now attempt to give an interpretation of the re- 
duced cluster dimension in terms of intragalactic gas 
porosity and star formation rate. The volume filling fac- 
tor / of the empty or low-density regions of a self-similar 
medium, the porosity, can be related to the medium's frac- 
tal dimension in three dimensions, D3, by 



/ 



1-1^ 



3--D3 



where r; and r„ arc the lower and upper boundary of the 
relevant scaling range (e.g., Turcotte 1992). From Fig. 7, 
left, and as mentioned above, we infer r; « 100 pc and 
Tu ~ 1000 pc. This approach to galaxy porosity is anal- 
ogous to Elmegreen's (1997) treatment of fractal inter- 
stellar gas clouds, the porosity of which was character- 
ized by fiCM = 1 — C*^^^/"^-*^^, with a maximum density 
contrast of C ~ lO'^ — 10^ for the intracloud gas. The 
two approaches are formally and numerically similar if we 
identify C = (ru/ri)^ w 10"^. Qualitatively, dwarf irregu- 
lar galaxies may thus be considered as huge star-forming 
clouds similar to fractal intragalactic star- forming clouds. 
Solving for the dimension, we obtain 



D3 « 3 + log{l-f). 



(3) 



Interpreting Fig. 7, right, in terms of galaxy porosity, we 
have to take into account that the scaling dimension of 
a projected isotropic sclf-similiar object is one less than 
the true dimension (Elmcgrcen & Elmegreen 2001), thus 
D3 = D + 1. Wc then learn that on average fainter 
galaxies with on average lower cluster dimensions, i.e. 
with stronger clustering properties, are also more porous 
(D « 1.5, / « 0.7) than brighter galaxies {D « 1.9, 
/«0.2). 

Theoretically porosity is thought to be crucial for the 
self-regulation of disks, and one expects an increasing 
star-formation rate to be accompagnied with decreas- 
ing porosity (Silk 1997, cqu. 7). This holds empiri- 
cally as well, as we will sketch now. For dwarf irregu- 
lar galaxies the area-normalized star formation rate is 
correlated with the galaxy's extrapolated central sur- 
face brightness: from Fig. 7a in van Zee (2001) we in- 
fer jjPg « -1.79±o.i8 \og{SFR/area) + 18.214±o.334, with 
area = 7r(1.5i?d)^ and Rd being the exponential-model 
scale length in kpc. On the other hand, an ordinary least- 
squares bisector fit (Isobe et al. 1990) to the data of Fig. 7, 
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Fig. 8. Comparison of three morphological indices, as ap- 
plied to bright lumps within irregular dwarf galaxies: the 
lumpiness index X, the normalized concentration index 
CI{x) (where x = B.25/Rd), and the reduced cluster di- 
mension D (taking into account the effect of differing disk 
scale lengths). 



right, yields /ig = — 3.25±i.oi D + 26.862-ti.5(;o. shown as 
line in the figure. Equating the two expressions, inserting 
equation (3), and remembering D = D3 — 1, we finally 
deduce 



SFR [Mgyr-i] « 0.45 (1 - f)'-^ [Ra [kpc])^ 



(4) 



Within our model treatment of dwarf irregular galaxies 

being self-similar objects we thxis have semi-cmpirically es- 
tablished a statistical relation between SFR, scale length, 
and porosity, in the sense that for a given scale length 
galaxies with higher SFRs arc also less porous. Note that 
for a given scale length, equation (4) predicts a maximum 
SFR. However, porosity as defined above has to be un- 
derstood as a conceptual parameter and not as a quantity 
describing reality in detail. The parameters possibly influ- 
encing the mean porosity of a galaxy are manyfold (gas 
density, gas pressure or velocity dispersion, gas metallic- 
ity, supernova energy release), forming an intricate, inter- 
dependent parameter set (Silk 1997). 

6. Comparing clustering, spreading, and relative 
luminosity of the lumps 

In the course of this paper we encountered three indices, 
namely the lumpiness index X (Section 3), the concentra- 
tion index CI (Section 4.3), and the cluster dimension D 
(Section 5). The first one has been introduced only, but 
has not been applied until now. To make up leeway and 
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to draw some first conclusions, we plot in Fig.|g the val- 
ues resulting from applying the three indices against each 
other. Several distinctive features are present. First, and 
most notable, the lumpiness index X anticorrelates with 
cluster dimension (upper left panel): the higher the lumpi- 
ness index , i.e. the higher the fraction of total galaxy light 
in the lumps, the smaller the cluster dimension, i.e. the 
more clustered are the bright lumps. Without plotting the 
corresponding graphs, we note that the lumpiness index 
shows no dependence on galaxy luminosity; this is con- 
trary to the behaviour of the cluster dimension. In terms 
of porosity, this behaviour means that lumpier galaxies are 
also more porous, and vice versa. 

Second, the degree of lumpiness and concentration be- 
have such that the galaxies with the highest values for the 
lumpiness index all have low values for the concentration 
index (lower panel). In other words, high fractional lump 
luminosities come with star-forming complexes that are 
widely scattered within a galaxy disk. On the other hand, 
those galaxies with lumps very centrally concentrated, i.e. 
with C/>3, come with low fractional lump luminosities, 
i.e. with low X values. This is somewhat surprising because 
there are BCD- and BCD-like galaxies with central star- 
bursts (high CI) that are expected to be very lumpy (high 
X ) at the same time. But indeed, no galaxy with high 
fractional lump luminosity is observed, not even a BCD 
galaxy, to be centrally concentrated. This means that even 
for actively star-forming galaxies the main body of light 
is still clearly dominated by the total galaxy light. We dis- 
claim a selection bias to the disadvantage of dwarfs with 
central starbursts in the next section by providing two 
representative examples. 

Finally, no relation seems to hold between concentra- 
tion index and cluster dimension (right panel) . In particu- 
lar, galaxies with highest lump concentrations present any 
cluster dimensions. In other words, lump location and the 
degree of self-similar clustering are independent of each 
other; this argues, as done in Section 4, for the introduc- 
tion of some mean galaxy porosity that may vary among 
galaxies. 

7. Discussion and conclusions 

Regarding the azimuthally integrated, radial distribu- 
tion of bright lumps — corresponding to star-forming 
complexes — in dwarf irregular galaxies, we find them 
non-uniformly distributed. While in individual galaxies 
the number distribution is non-monotonic and rugged, 
the summed-up distribution for all galaxies of our sam- 
ple manifests the hidden constraint, which is a re~^- 
distribution closely tracing the underlying older popula- 
tion. More precisely, in terms of radial number density dis- 
tribution the lumps follow an exponential decay with scale 
length about 10 percent smaller on average than that of 
the blue continuum light. This is consistent with studies of 
the radial distribution of H II regions in intermediate- type 
spiral galaxies (Athanassoula et al. 1998). The fact that 
each component of the average disk — from star-forming 



site number density to surface brightness of the total light 
— is approximately exponential is a hint that luminous 
exponential disks are born rather than made, consistent 
with the accretion scenario for the viscous evolution of 
galaxy disks (e.g., Ferguson & Clarke 2001). Star-forming 
complexes in irregular dwarf galaxies can be found out to 
large radii, as already emphasized by Schulte-Ladbeck & 
Hopp (1998), Brosch et al. (1998), and Roye & Hunter 
(2000). The presence of a tail in the accumulated radial 
number distribution of star forming regions out to at least 
six optical scale lengths (Fig. 4) indicates that the dis- 
tributions of dwarf irregulars are truncated at rather low 
gas density thresholds for star formation. ( van Zee et al. 
1997, Hunter et al. 1998, Pisano et al. 2000); this seems 
to be different with many spiral galaxies for which sharp 
to weak truncations, starting at galactocentric radii of 2-4 
near-infrared or 3-5 optical scale lengths, have persistently 
been reported (e.g., recently, Florido et al. 2001, Kregel et 
al. 2002, Pohlen et al. 2002). 

Beside the presence of main or primary peaks in the 
radial lump number distribution at slightly less than one 
optical scale length on average, there is — contrary to 
the expectations for exponential-disk systems — the fre- 
quent occurence of secondary peaks at about two scale 
lengths. As simple simulations show, it is consistent with 
the idea of triggered star formation based on a stochas- 
tic self-regulation scenario. However, some of the brighter, 
larger galaxies exhibit pronounced primary peaks at two 
scale lengths and show minor, secondary peaks around one 
scale length. For these galaxies with a reversed peak pat- 
tern the simulations indicate that shear-induced star for- 
mation around the disk's turnover to differential rotation 
could be at work; we feel this issue worth a deeper investi- 
gation: given a lump statistics based on higher-resolution 
images and linked to detailed rotational velocity data, 
and possibly supplemented with information on large-scale 
magnetic field structures, this may lead to some subtle 
but decisive insights related to star formation in irregular 
galaxies. Also, the possible role of bars or bar-like central 
features should be carefully considered. However, as none 
of our four galaxies with principal peaks appering around 
two scale lengths is classified as "barred", and because 
Roye & Hunter (2000) did not see a preferential location 
of H II regions towards the ends of bars in the two can- 
didate galaxies of their sample, we do not consider this 
mechanism as being effective in shaping the number dis- 
tribution of lumps. 

The observation that the scale lengths are the larger 
the older the underlying respective population is, goes in 
hand with the finding that in star-forming dwarf galax- 
ies the oldest populations are also the most extended ones 
(Gardiner & Hatzidimitriou 1992, Minniti & Zijlstra 1996, 
Minniti et al. 1999, Harris & Zaritsky 1999). The com- 
mon interpretation is that of an age-related dispersion of 
stars. Because dynamical disk heating tends to saturate at 
a fixed velocity dispersion (Freeman & Bland-Hawthorn 
2002) the amount of radial spread introduced by dynami- 
cal heating is expected to be larger in smaller galaxies with 
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lower gravitational binding energies (J. Gallagher, private 
communication). Indeed, the fact noted in Parodi et al. 
(2002) that with increasing scale length the disk color gra- 
dients become systematically less positive and even start 
being weakly negative when going from dwarf irregular 
to low-surface brightness and spiral galaxies — which is 
equivalent to larger galaxies having red-to-blue band scale 
length ratios below one — seems to support the idea that 
the extent of the red component decreases as galaxy mass 
increases. However, for spirals part of this colour trend is 
most certainly a metallicity effect. 

We provide no investigation of the azimuthal lump 
distribution; where necessary we assumed axisymmetry. 
Applying different asymmetry indices to dwarf and nor- 
mal irregulars. Heller et al. (2000) and Roye & Hunter 
(2000) obtained opposing results that were dependent on 
the definition of the chosen index (asymmetry as a ques- 
tion of perception). Also, it was found that while the con- 
veniently applied rotational asymmetry index may be used 
as a first discriminator between (distant) elliptical and spi- 
ral/irregular galaxies (Schade et al. 1995, Conselice et al. 

2000) , this parameter actually is strongly dependent on 
recent star formation (Takamiya 1999, Mayya & Romano 

2001) and thus must be correlated with the lumpiness in- 
dex as adopted in the present paper. As we have no data 
for the star formation rate of our galaxies, an immediate 
comparision cannot be offered, however. 

Applying a concentration index CI that is normal- 
ized according to the exponential-disk structure of a mean 
lump distribution leads to consistent results for varying 
aperture sizes. It may also remove the discrepancy found 
by Heller et al. (2000) between the CIho values for actual 
galaxies and the lower ones for simulated galaxies with 
random star formation region positions. Roye & Hunter 
(2000) pointed out an increased scatter of concentration 
indices for faster rotating galaxies of their sample; we no 
longer see this effect with our larger sample. 

Comparing concentration {CI) with lumpiness (X ) we 
find the galaxies with a high percentage (> 10%) of light 
stemming from the lumps showing low to moderate con- 
centrations {CI 5 2), i.e., galaxies with lumps that are 
widely scatterd within the disk maintain a higher fraction 
of the total B luminosity (Fig. ^). On the other hand, 
for very concentrated galaxies {CI ^ 2) less than about 
10 percent of the light is due to the lumps; actually, for 
these galaxies the values for the lumpiness index X fluc- 
tuate around the sharp value of 7% attributed to most 
of the (barred) spiral and irregular galaxies observed by 
Elmegreen & Salzer (1999). They suggested the similar- 
ity of the blue-band light fraction in complexes for several 
galaxies of different Hubble types and different total lu- 
minosities being due to similar star formation efhciencies. 

It remains, however, remarkable that galaxies with 
very high lump concentrations arc not among the galaxies 
showing high B luminosity fractions of the lumps (Fig. 
8, bottom panel). One may wonder whether some of the 
nearby, well-known BCDs with central starbursts would 
agree with this conclusion as well, i.e. whether their cen- 



tral bursts should not dominate the total light content. 
We therefore examined where the typical starburst dwarf 
irregular galaxies NGC 1569 and NGC 1705 would fill in 
the CI vs. X diagram. Based on the 48 brightest central 
star clusters of NGC 1569 as compiled in Hunter et al. 
(2000) and adopting a distance of 2.5 Mpc and a galaxy 
absolute magnitude of My = —17.99 mag, we estimate a 
mere X < 0.10 for NGC 1569. Similarly, for the super-star 
cluster dominated amorphous galaxy NGC 1705, the abso- 
lute magnitude for the 7 brightest clusters is about My = 
— 14.1 according to the data in O'Connell et al. (1994), 
whereas the galaxy has a magnitude of My = —16.13 at a 
distance of 5 Mpc, implying only X ~ 0.15. Thus indepen- 
dent of the exact concentration indices for the brightest 
clusters, these two small but representative galaxies with 
central starbursts would not occupy the empty part of the 
C/ vs. X diagram. This means that the empty region in 
this figure is not a selection effect, but may be related to 
our adopted procedures in determining the corresponding 
indices. 

While the concentration index is a measure for lump 
spreading, the cluster or correlation dimension provides 
information on the scaling behaviour for lump-to-lump 
distances. We found the lump cluster dimensions — cor- 
rected for the effect of radial abundances according to the 
annulus-integrated exponential distribution — to lie be- 
tween 1.3 and 2.0 and to gently correlate with extrap- 
olated central surface brightness and absolute magnitude 
of the host galaxy. At the same time the cluster dimension 
is weakly anticorrelated with the lumpiness index. 

Cluster dimension (or porosity) as introduced in this 
paper may be intimately linked to the sizes of the largest, 
kiloparsec-sized lump compounds (Elmegreen et al. 1996) 
or to the sizes of star-forming, collapsed expanding shells 
(e.g., Walter 1999). Elmegreen et al. (1996) found the 
sizes of the largest compounds within spiral and irregu- 
lar galaxies to approximately scale with the square root 
of the galaxy luminosity, or, if normalized by the galaxy 
semi- major axis R25 , small galaxies to have slightly larger 
relative compound diameters than larger galaxies. This 
was hypothesized to result from gravitational instabilities 
with the Jeans length or the mean virial density varying 
with galaxy luminosity. Walter (1999) suggests the holes 
in dwarf irregulars to be larger than those in late-type 
spirals because small galaxies have lower masses and cor- 
respondingly lower gravitational potentials and lower am- 
bient ISM gas densities, favoring HI shells to grow larger. 
Alternatively, we are temptatively interpreting the clus- 
ter dimension of a galaxy in terms of the volume filling 
factor of empty regions in a fractal medium, i.e. in terms 
of porosity (defined analogous to Elmegreen 1997), and 
find the following statistical trends: (i) fainter galaxies 
tend to be more porous; (ii) more porous galaxies also have 
a lumpier morphology with lower central lump concentra- 
tions; (Hi) the more porous the galaxy, the lower the star 
formation rate per kp(? (equation 4-)- While these trends 
are not unexpected, we provide an objective and quantita- 
tive statistical treatment of these. 
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Fig. 9. SSPSF simulation of an exponentiated disk with 80 corotating rings and a transition from rigid to strong 
differential rotation at two scale lengths (ring 32, solid line), after 500 time steps. The three panels show the pattern of 
newly activated cells (filled circles) together with those created less than 10 time steps ago (dots), the radial number 
distribution of the active cells (bin width is 8 rings or half a scale length), and the corresponding number density 
distribution. Note that for the number distribution the peak at one scale length (ring 16) — the presence of which is 
expected for non- differentially rotating exponential disks — is markedly surpassed by a peak around two scale lengths. 



Porosity, or self-similarity, as observed with the 
bright-lump distribution within dwarf irregular galaxies 
reflects the sclf-rcg\ilatcd cvohition of the interstellar 
medium, with stellar feedback and self-gravitation being 
the main mechanisms. It is thus not to be confused with 
the still self-similar pattern of dispersed stellar aggregates 
that initially formed from the fractal interstellar gas, 
obeying the canonical value of D=1.3 for turbulence- 
driven star formation (Elmegreen & Elmegreen 2001). 
Our sub-galactic areas are showing D>1.7, i.e. much 
higher fractal dimensions. We suspect that in dwarf 
irregular galaxies either the dispersive redistribution of 
stars is indeed much more effective than in spiral galaxies 
(Elmegreen & Hunter 2000), and/or that feed-back 
regulation is responsible for the partially randomized 
position patterns (e.g., due to the intersection of giant 
shells). 
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Appendix A: The SSPSF model 

The implementation of the stochastic self-propagating star 
formation (SSPSF) disk model is a variant of the prescrip- 
tions found in the reviews of Schulman & Seiden (1986) 
and, in more detail, Seiden & Schulman (1990). Major 
modifications are an exponentiated initial disk structure, 
a transition from a linearly increasing to a flat rotation 
curve, and the inclusion of spontaneous star formation. In 
particular, our model consists of 80 corotating rings, each 



with 6R cells, R being the ring number or radius. Thus 
we have a total of 18960 cells. Initially, for each ring the 
number of occupied cells, N{r), has a probability propor- 
tional to a Gamma distribution Rexp{—R/Rd), with the 
scale length taken to be Rd = 16. On average the galaxies 
of our sample have a scale length of about 0.7 kpc, thus 
the linear size of each cell corresponds to about 44 pc. 
The polar angles for the N{r) cells in a ring are randomly 
chosen. At the start there arc a total of 300 occupied cells. 
The life time of an occupied cell is 10 time steps, corre- 
sponding to about 10 Myr. The first time step following 
the activation of a cell, an empty neighbouring cell has a 
probability of 0.21 to become activated too. This imple- 
ments the idea that the stellar wind of massive stars cre- 
ated in a cluster travels with a wind velocity of about 40 
kms~^ for about 1 Myr forming an increasingly dense shell 
that eventually fragments or hits other overdense regions, 
thus giving chance to the formation of new star-forming 
sites. At each time step, an additional 20 new cells are 
spontaneously activated; this sustains the number-density 
profile being exponential. Typically, an equilibrium occu- 
pation of around 1400 cells is reached after a few dozen 
time steps (providing a filling factor around 0.07), 10 per- 
cent of which having been just activated, and with the 
number density distribution remaining rather exponential 
out to about four scale lengths; further out there is a rapid 
drop or tnmcation of star- forming regions. 
Either the rings rotate rigidly, with circular velocity be- 
ing proportional to the ring radius, or differential rotation 
may be imposed by additionally demanding a constant cir- 
cular velocity beyond a turnover radius Rt = 2i?d. For the 
simulations that included shear a flat rotation curve ve- 
locity of 500 kms"'^ was enforced, which is an unrealistic 
factor of ~10 faster than observed for a typical dwarf ir- 
regular. However, this was only to clearly demonstrate the 
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effect of strong shear within the model. As long as the cir- 
cular velocity is much higher than the radial propagation 
velocity, shear seems to increase the star formation rate. 
For each simulated galaxy the number (and number den- 
sity) distribution after 500 time steps was stored for use 
in the study of Section 4.4. The lump patterns and their 
corresponding distributions for a typical simulated galaxy 
can be seen in Fig. 9. Note that we simulate only the oc- 
curence of new and the presence of young stellar clusters, 
but otherwise assume a prevailing underlying population 
of older stars. 
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